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ABSTRACT 

Strong gravitational lensing has traditionally been one of the few phenomena 
said to oppose a large cosmological constant; many analyses of lens statistics have 
given upper limits on that are marginally inconsistent with the concordance 
cosmology. Those conclusions were based on models where the predicted number 
counts of galaxies at moderate redshifts (z ~ 0.5-1) increased significantly with 
f^A- I argue that the models should now be calibrated by counts of distant 
galaxies. When this is done lens statistics lose most of their sensitivity to the 
cosmological constant. 

Subject headings: cosmological parameters — gravitational lensing — galaxies: 
evolution 



1. Introduction 

Popular opinion seems to have settled on a "concordance" cosmology dominated by 
dark energy. The cosmic microwave background indicates a flat geometry (e.g., de Bernardis 
et al. 2000; Hanany et al. 2000; Pryke et al. 2002), cluster mass-to-light ratios indicate a 
low matter density (e.g., Carlberg, Yee & Ellingson 1997; Bahcall et al. 2000), and type 
la supernovae indicate cosmic acceleration (e.g., Riess et al. 1998; Perlmutter et al. 1999). 
The popular cosmology has matter content Qm ~ 0.30-0.35 and dark energy Qx ~ 0.65- 
0.70 such that Qm + &x = 1- For simplicity I assume that the dark energy corresponds to 
a cosmological constant A, although my analysis could easily be extended to quintessence 
models (e.g., Waga & Miceli 1999). 
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One phenomenon that has traditionally stood out from the concordance model is strong 
gravitational lensing. The statistics of strong lenses are sensitive to the cosmological param- 
eters via the cosmological volume element (e.g., Turner 1990; Fukugita, Futamase & Kasai 
1990), and analyses of the data have yielded upper limits on the dark energy at the level of 
Qa < 0.66 at 95% confidence (e.g., Kochanek 1996; Falco, Kochanek & Muhoz 1998). There 
are small systematic uncertainties in the upper limit due to assumptions about the lens 
sample and the amount of dust extinction in lens galaxies (e.g., Falco et al. 1998; Helbig et 
al. 1999; Waga & Miceli 1999; Cooray, Quashnock & Miller 1999). Larger systematic effects 
arise from uncertainties in the local luminosity function of galaxies. Chiba & Yoshii (1999) 
argue that by adopting luminosity functions from different surveys they can relax the upper 
limit on Q\ and even find models that favor values in the range Q\ ~ 0.5-0.8. Kochanek et 
al. (1998) respond by acknowledging the systematic uncertainties but defending their choice 
of luminosity function as the one that is most consistent with the observed luminosities of 
lens galaxies. This controversy will soon be resolved by new measurements of the luminosity 
function from the SDSS and 2dF surveys that appear to eliminate most of the traditional 
problems (Blanton et al. 2001; Cross et al. 2001). 

While debating the details, the previous studies agreed on the idea that raising Q\ 
dramatically increases the expected number of lenses; they agreed on the trend and only 
contested the zero point. In this Letter I question the trend itself, based on number counts 
of distant galaxies. I argue that the lens statistics models used to obtain upper limits on 
f^A are inconsistent with galaxy counts at z ~ 0.5 — not with any particular value of the 
counts, but with the general idea that they can be measured and used as constraints on the 
models. I modify the models to be calibrated by galaxy counts and consider how the new 
models depend on Q\. For simplicity I consider only flat cosmologies. 



where dn/dM is the mass function of deflectors, A(M, zi, z s ) is the cross section for lensing 
by a deflector of mass M at redshift zi, and D(z) is the comoving distance. When analyzing 
real lens data the optical depth must be modified to account for magnification bias, or the 
fact that lensing magnification effectively changes a survey's flux limit (e.g., Turner, Ostriker 
& Gott 1984). Magnification bias is largely insensitive to cosmological parameters, so for 
our purposes the simple optical depth is adequate. 



2. 



Fundamentals 



The optical depth for lensing of sources at redshift z s is 




(1) 
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Nearly all arcsecond-scale lenses are produced by detectable galaxies, so many analyses 
of lens statistics model the deflector population with observed galaxy populations. The 
luminosity function of galaxies is often parameterized as a Schechter function with the form 
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where is a characteristic luminosity, is a characteristic comoving number density, and a 
is the power law slope at the faint end. Luminosity must be converted to mass for the lensing 
analyses, and this is usually done with the empirical Faber-Jackson relation L/L* = (cr/o"*) 7 
where a is the velocity dispersion. 

Luminosity naturally changes with redshift due to passive evolution: stellar populations 
fade with time as stars die out. The number density may change with time due to mergers. 
The number evolution can be parameterized in various ways, but Lin et al. (1999) propose 
the form 

n t {z) =n,(0)10°- 4P2 , (3) 

so that the number evolution parameter P can be combined with another parameter rep- 
resenting luminosity evolution expressed in observers' units of magnitudes. I adopt this 
form since the exact parameterization is not very important for my analysis. Other than 
luminosity and number evolution, we may assume for simplicity that a, 7, and cr* are all 
independent of redshift (see Lin et al. 1999). 

Lens galaxies are usually approximated as isothermal spheres, which is consistent with 
lens data, galaxy dynamics, and X-ray ellipticals (e.g., Fabbiano 1989; Rix et al. 1997; Treu 
& Koopmans 2002). The lensing optical depth can then be written as (e.g., Kochanek 1993) 
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where D oh D os , and D ls are comoving distances between the observer, lens, and source, 
written here in units of the Hubble distance r# = c/H . 



3. Results 

If there is no number evolution (P = 0), then eq. (4) yields the standard result for the 
lensing optical depth (e.g., Kochanek 1993), 

t(z.) = ^ Dl , (6) 
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where the optical depth is proportional to the cosmological volume D^ s and hence very 
sensitive to as shown in Fig. 1. If we allow for some number evolution but assume 
that n*(z) is identical for all cosmologies (as done by Jain et al. 2000), the optical depth 
curves shift slightly but retain the same qualitative behavior. In other words, if we assume 
the same number density (whether constant or not) for all cosmologies then we recover the 
standard result that the number of lenses increases significantly as Qa increases, with little 
sensitivity to the actual amount of number evolution. The increase is large enough that 
previous analyses of lens statistics were able to obtain interesting upper limits on £1\ even 
from small lens samples. 

Although the focus of the models is the lensing optical depth, a corollary prediction is 
the number of galaxies per unit redshift per unit area on the sky, 2 

- = »,WD( 2 f-r(i + «). (7) 

Fig. 2 shows the predicted number counts for the models from Fig. 1. If n*(z) is the same 
in all cosmologies then number counts of galaxies at moderate or high redshifts are very 
sensitive to Q\. For example, at z — 0.5 models with Q\ = and 0.7 differ by a factor 
of two. This is the basis of the classic idea to use galaxy number counts to constrain the 
cosmology (see Sandage 1997; Driver 2002). Lens statistics boil down in some sense to 
galaxy counts, where the objects being counted are lens galaxies (mostly massive ellipticals 
at z ~ 0.3-1). 

Galaxy counts at moderate redshifts (z ~ 0.5-1) can now be measured directly (e.g., 
Lin et al. 1999; Davis et al. 2001; Eisenstein et al. 2001), so we can imagine using them to 
test the models. Let us consider not any particular value for the counts, but just the general 
idea that they can be measured. If the counts at z ~ 0.5-1 are measured then at most one of 
the curves in Fig. 2 can be correct. What happens if instead we insist that all models agree 
with some observed galaxy counts dN/dzl In an idealized situation where number counts 
could be measured well at all redshifts, the optical depth would be 

^ ' \c) r(l + a) J \D os J dz 1 ; 

In other words, the optical depth would simply be an integral over the observed number 
counts, weighted by the distance ratio D ls /D os . The direct dependence on the cosmological 
volume D% s would disappear, leaving only a weak dependence on cosmology through the 
distance ratio. 



2 For simplicity this analysis omits luminosity cuts, but it would be straightforward to add them. 



- 5 - 



In practice the number counts will be measured well for some finite range of redshifts. 
Given counts dN/dz at redshift z ohs and parameterizing the evolution as in eq. (3), we would 
infer a number evolution parameter 3 



Note that this value depends on the assumed cosmology through the volume element D 2 dD/dz. 
Assuming two different cosmologies ("1" and "2") would lead to two different values of P 
that are related by 



Thus, when models are calibrated by observed galaxy counts at moderate redshifts, the 
function n*(z) cannot be the same for all cosmologies. The number density is the number 
counts divided by the cosmo logical volume element, so it necessarily depends on the assumed 
cosmological parameters. By using the same n*(z) in all cosmologies, traditional models for 
lens statistics have violated this condition. 

Let us construct new models that do not violate it. Imagine that from observed galaxy 
counts we infer a number evolution parameter P in a cosmology with ^a = 0. For other 
values of Qa we can use eq. (10) to compute the self-consistent value of P, as in Fig. 3. 
Models with different values of Q\ now agree on the predicted number counts at z ^ s (by 
construction) and are very similar over a wide range of redshifts z < z ohs , as shown in Fig. 4. 4 
In other words, the new models are consistent in that they agree on an observable quantity, 
galaxy number counts, independent of 

Fig. 5 shows the lensing optical depth as a function of Qa for these new models. The 
difference from the old models is striking: the dependence of the optical depth on Q\ is 
much smaller than before. In the old models, going from Q\ = to Q\ = 0.7 increased 
the number of lenses by a factor of 2.7-3.3 (depending on the value of P), while in the new 
models the increase is only 40-43% if the calibration is at z ^, s = 1.0, or a mere 5-16% if 
z bs = 0.5. Models constrained by galaxy counts at all redshifts as in eq. (8) (not shown) lie 
between the z ^ s = 0.5 and z ^, s = 1.0 curves. 

The reduced dependence on is a direct result of calibrating the models with observed 
counts of distant galaxies. This qualitative conclusion does not depend on any particular 



3 If counts were known over a range of redshifts, this expression and eq. (10) would include appropriate 
averages over z bs- 

4 The various models could be made to match over a wider range of redshift by introducing more param- 
eters to describe the evolution. 






(10) 



- 6- 



values I have assumed, but arises simply from the idea that models for lens statistics can 
and should be constrained to agree with measurable galaxy counts at redshifts relevant for 
lensing. There is some residual sensitivity to cosmology through the distance ratio Di s /D os 
(see eq. 8), but it is much smaller than the volume effect and thus will be harder to detect. 

The prospects for using real galaxy counts to calibrate models for lens statistics are good. 
The CNOC2 field galaxy redshift survey already contains ~5000 galaxies at 0.12 < z < 0.55 
(Lin et al. 1999), the SDSS Luminous Red Galaxy sample will have redshifts for ~100,000 
early-type galaxies out to z ~ 0.5 (Eisenstein et al. 2001), and the DEEP redshift survey 
will include ~60,000 galaxies at z > 0.7 (Davis et al. 2001). The calibration of the models 
will be limited to some extent by completeness and selection effects, but with large samples 
it will be possible to make different cuts on the data to understand those effects (see Blanton 
et al. 2001). The calibration will also be limited by our understanding of the relationship 
between luminosity and mass — but that has always been true for models of lens statistics. 
The best bet is probably to turn the problem around: rather than trying to constrain f^A, 
a joint study of galaxy counts and lensing could attempt to distinguish between evolution 
in luminosity and evolution in mass. That would be perhaps the most interesting future 
application of lens statistics. 



4. Conclusions 

The optical depth for lensing is basically proportional to the number of galaxies on the 
sky at redshifts z\ ~ 0.3-1. Traditionally that number was not well known, so models for lens 
statistics adopted a number density of galaxies and multiplied by the cosmological volume 
to get the number. With the assumed number density held fixed, the expected number of 
lenses was proportional to the cosmological volume and hence very sensitive to Q\. 

Number counts of distant galaxies can now be measured directly, and they are inconsis- 
tent with the idea that the number density is independent of Q\. This general point holds 
whether there is much or little redshift evolution in the galaxy population. Constraining 
models for lens statistics to agree on distant galaxy counts makes them far less sensitive 
to f2 A . Whereas the old models saw a factor of 3 change in the number of lenses between 
f2 A = and f2 A = 0.7, in the new models the change is <10-40%. Using lens statistics to 
constrain f2 A m &y still be possible, but will be difficult. 

This paper has focused on models where the deflector population is derived from ob- 
served galaxy populations, which I refer to as phenomenology models. In an alternate class 
that I call theory models, the deflector population is described with a mass function from 
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structure formation theory; the resulting models are sensitive to cosmology not only through 
the volume element but also through Qm and the growth of structure. In theory models 
the predicted number of lenses decreases as increases (for flat cosmologies; Porciani & 
Madau 2000; Li & Ostriker 2002), which strongly disagrees with old phenomenology models 
but is less different from my new models. In quintessence models, phenomenology and the- 
ory models do agree that making the equation of state more negative increases the predicted 
number of lenses (Waga & Miceli 1999; Sarbu, Rusin & Ma 2001). The differences between 
phenomenology and theory models clearly need further study. Nevertheless, I would argue 
that the focus of lens statistics should move away from constraining JIa and toward learning 
about the population of dark matter halos out to z ~ 1. 

I am very grateful to the referee, David Rusin, for important criticisms and insightful 
comments. I also thank Huan Lin for interesting discussions and Wayne Hu for comments 
on the manuscript. This work was supported by NASA through Hubble Fellowship grant 
HST-HF-01141.01-A from the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5- 
26555. 
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Fig. 1. — The lensing optical depth as a function of Oa> for a source redshift z s = 2. The 
optical depth is normalized by its value at f^A = 0. The number evolution parameter P is 
assumed to be the same for all cosmologies. 
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Fig. 2. — Galaxy number counts for the P = models from Fig. 1. Different line types 
indicate different values of Q\. 
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Fig. 3. — The change in the number evolution parameter P with if the models are 
calibrated by galaxy number counts at redshift z ^, s = 0.5 (solid) or ^ b s = 1.0 (dotted); see 
eq. (10). The curves correspond to P — 1 (top), (middle), and —1 (bottom). 
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Fig. 4. — Galaxy number counts for new models calibrated by some observed counts at 
z obs = 0.5 (top) or z ohs = 1.0 (bottom). Different line types indicate different values of Qa- 
The hypothetical measured counts correspond to P — 1, 0, and —1 for the top, middle, and 
bottom set of curves in each panel. Note that in the limit z Q \, s — > we would recover Fig. 2. 
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Fig. 5. — Similar to Fig. 1, but with results for new models calibrated by galaxy counts at 
redshift z Q h s = 0.5 or 1.0. 
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